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Abstract 

An analytical method based on the classical ruin problem is developed to compute the mean 
reaction time between two walkers undergoing a generalized random walk on a Id lattice. At each 
time step, either both walkers diffuse simultaneously with probability p (synchronous event) or one 
of them diffuses while the other remains immobile with complementary probability (asynchronous 
event). Reaction takes place through same site occupation or position exchange. We study the 
influence of the degree of synchronicity p of the walkers and the lattice size N on the global 
reaction's efficiency For odd N, the purely synchronous case (p = 1) is always the most effective 
one, while for even N, the encounter time is minimized by a combination of synchronous and 
asynchronous events. This new parity effect is fully confirmed by Monte Carlo simulations on Id 
lattices as well as for 2d and 3d lattices. In contrast, the Id continuum approximation valid for 
sufficiently large lattices predicts a monotonic increase of the efficiency as a function of p. The 
relevance of the model for several research areas is briefly discussed. 
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1. INTRODUCTION 



Recently, the modelling of diffusion-controlled reactions on lattices has attracted renewed 
interest due to synergies with some research areas like e.g. heterogeneous catalysis [H 0|, 
trapping problems 0], spin models [H, game theory [H|], population dynamics 6] or biological 
problems 0]. An issue common to these systems is the important role played by the coex- 
istence of different intrinsic time scales, the lattice characteristics (size and dimensionality) 
and many-body effects. The interplay of these ingredients may strongly affect the efficiency 
with which statistical processes such as front propagation on catalytic substrates, the spread 
of an infection, kink propagation in magnetic systems, or exciton trapping in photosynthetic 
cells take place. In order to obtain analytical insight for such systems, it is often necessary 
to make simplifying assumptions which nevertheless preserve their main generic features. 

In this spirit, a prototypical diffusion- reaction system consists of two interacting walkers 
performing nearest neighbor jumps on a lattice The walkers are assumed to react with 
each other whenever they meet at the same lattice site or attempt to exchange positions. 
Each of such "encounters" results in instantaneous reaction, the process therefore being 
diffusion-controlled. Regardless of the particular outcome of the reaction, a measure for its 
efficiency will be clearly given by the mean encounter time of both walkers. 

In a pioneering work, Montroll investigated a simplified version of the problem in which 
one of the walkers is stationary, thereby playing the role of a fixed trap [9]. In a recent 
article Montroll's results have been extended with the help of a Markovian method to 
account for the possibility of simultaneous displacement of the walkers. In the present work, 
the problem is further extended to the case in which the motion of the walkers consists of 
a random sequence of two different events: at each tick of the clock, a synchronous event 
takes place with probability p, i.e., both walkers hop simultaneously to randomly chosen 
nearest neighbor sites. Alternatively, an asynchronous event takes place with probability 
1 — p, i.e. one of the walkers performs a nearest neighbor jump while the other remains 
immobile. Thus, the parameter p interpolates between the one walker plus trap case studied 
by Montroll (p = 0) and the case of two simultaneously moving walkers (p = 1). A central 
question we shall investigate here concerns the influence of the degree of synchronicity p of 
the walkers and the size N of the lattice in which they are embedded upon the reaction's 
efficiency. 

Our model may be of interest in several contexts. One can e.g. easily accomodate it to 
allow for events in which both walkers remain immobile in the original reference frame. Thus, 
three qualitative different joint events become possible for the two-walker system, which can 
e.g. be interpreted as resulting from a combination of two internal states for each of the 
walkers, namely a diffusing and an immobile state. Such two-state random walks 13, 3 
are frequently used to model chromatographic [14] or electrophoretic separation processes, in 
which the propagation of charged particles in an external field may be occasionally stopped 
by entanglement with the molecules of the substrate. Such systems, along with hopping 
models for transport and recombination of carriers in solids (Til lla 17], might provide 
additional motivation for considering a generalized version of the random walk. 

A different approach suggests to regard the fluctuations in the diffusivity of the walkers 
as a result of random fluctuations of lattice sites switching between a conducting and a 
stopping state, whereby the translational invariance of the lattice is preserved on average. 
An alternative formulation of the problem in terms of fluctuating dichotomous barriers 
between sites also seems possible [l8j . Such models for dynamic random media are relevant 
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for the description of several physiochemical processes like e.g. ligand diffusion in proteins 
[l9l | or proton migration in water jiol ]. 

On the other hand, if one goes back to the original picture of two dissimilar walkers A and 
B, the model may be regarded as a schematized starting point for describing the dynamics 
of exciton absorption in biological light-harvesting systems 0). In photosynthetic cells, a 
photon is absorbed by the pigment molecules (e.g. chorophyll) in the cell and may give 
rise to an excited energy state. The exciton hops by resonant energy transfer through 
a network or lattice of 200-500 pigment molecules (antenna system) and can be eventually 
trapped at a reaction center ji^], which is usually considered to be immobile within the time 
scale for trapping (a few hundred picoseconds). The exciton energy is then used to trigger 
a series of redox processes in the chain of chemical reactions leading to the production 
of sugars and carbohydrates. Thus, one of the time-limiting steps for the production of 
oxygen is precisely the absorption of the exciton by the trapping center. If we allow for a 
certain mobility of the reaction center (thereby generalizing Montroll's approach for exciton 
trapping), we can identify the latter with a slowly hopping walker, say the A walker, while 
the propagating exciton would play the role of the B walker. As long as the hopping 
rate of the A walker remains small, the situation may be identified with the almost purely 
asynchronous case (small p); this is normally the case in in vivo light-harvesting systems. 
However, a modification of the physical properties of the antenna system so as to make 
exciton propagation slower might, at least in principle, lead to interesting antiresonance 
phenomena as observed in our model. 

The paper is organized as follows. In section 2, we define the two- walker system in detail 
and show how it can be recast into an equivalent one-walker system with absorbing sites. 
In section 3, we report analytic work on the Id case and identify the principal differences 
between the purely asynchronous (p = 0), the purely synchronous case (p = 1) and the 
mixed case (0 < p < 1). These results are also supported by Monte Carlo simulations. In 
Section 4, similar results are found for the 2d and the 3d cases by means of simulations In 
section 5, the results are compared with the predictions of the continuum approximation 
valid for large lattices. Finally, section 6 summarizes the main conclusions and discusses 
possible extensions. 



2. FORMULATION OF THE PROBLEM: TWO- WALKER VS. ONE- WALKER 
PICTURE 

The starting point is a Id periodic lattice with N sites and discrete time dynamics (cf 
Fig. ^,). We place two walkers A and B on two distinct lattice sites and let them evolve at 
each time step as follows: 

1. with probability p, both walkers hop simultaneously to randomly chosen nearest neigh- 
bor sites (synchronous event). 

2. with probability 1— p, one of the walkers (no matter which one) remains at rest while 
the other performs a jump to a nearest neighbor site (asynchronous event). 

The walkers are assumed to be unbiased, i.e., their jumps are symmetric. We additionally 
assume that their jump directions are completely uncorrelated. An encounter takes place 
when both walkers 'land' on the same site or attempt to exchange positions. Each encounter 
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triggers instantaneously an irreversible reactive interaction, say the annihilation reaction 
A + B — > 0. The encounter time can thus be regarded as the characteristic reaction time 
governing the diffusion-controlled two-particle annihilation. 

In principle, the mean encounter time can be computed for a single initial configuration of 
the walkers; in many practical situations, though, one has little knowledge about the initial 
conditions. We shall therefore give preference to a definition of the mean encounter time 
which contains an additional coarse graining over all possible initial configurations. In the 
sequel, we shall denote this quantity by (n). The smaller (n), the more efficient the reaction 
will be. 

Note that, due to the translational invariance of the lattice, the physical distinguishability 
of the walkers is irrelevant for the computation of (n): in other words, it does not matter 
which of the walkers A or B hops more often, since (n) depends only on the relative motion 
of both walkers, the latter being fully characterized by p. Therefore, we shall assume for 
simplicity and without loss of generality that A and B are physically identical walkers and 
drop the labels A and B, as has been done in Fig. 

For p = 0, only reactions through simultaneous site occupancy are possible. As far as (n) 
is concerned, this case is equivalent to the one treated by Montroll [9|. However, if p > 0, 
reaction through position exchange becomes possible. This technical obstacle makes the 
analytical treatment of the problem more difficult. 

A first step to circumvent this difficulty is to take advantage of the translational invariance 
of the lattice and reformulate the problem in terms of a single walker. To do so, we must 
switch over to a new comoving reference frame in which one of the walkers remains stationary, 
thereby playing the role of a fixed perfect trap T (Fig. \2p). Clearly, an asynchronous event in 
the original two-walker system is equivalent to nearest neighbor hopping in the single-walker 
system, while a synchronous event results either in a walker's jump by two lattice sites (if 
both walkers hop in opposite directions in the original reference frame) or its remaining at 
rest (if they hop in the same direction). In this picture, reaction takes place any time the 
walker reaches or overreaches the trap T. 

We can now take advantage of the simple lattice geometry and the fact that the trap is 
perfect to unfold the iV-site lattice into an equivalent one with N + 1 sites and two perfect 
traps T sitting at each end site, as shown in Fig. |2]for N = 7. This transformation does not 
of course affect the characteristics of the walk; a walker jump in anticlockwise direction will 
be equivalent to a jump to the left by the same number of sites in the transformed lattice. 
Next, we replace each trap T by two (fictitious) reactive sites r (Fig. EJ). The random 
walk will instantaneously terminate when the walker lands on any of these r-sites. We shall 
therefore occasionally call the r-sites "absorbing" in what follows. 

Let us assign the coordinates 1 to iV — 1 to each of the non-absorbing sites, as shown in 
Fig. El Additionally, we term the r-sites at the left end —1 and and those at the right end 
N and N + 1. An attempt to overcome, say the left trap, triggered by a synchronous event 
will be equivalent to a jump from site 1 to site —1. In contrast, the walker's landing on the 
trap can be realized either by an asynchronous event (jump from site 1 to site 0) or by a 
synchronous event (jump from site 2 to site 0). If the dynamics is purely asynchronous, only 
one r-site at each end will be needed, since jumps by two sites are not possible in this case. 
Thus, events involving position exchange can be dealt with by going over to a one-walker 
picture and introducing fictitious absorbing sites. The mean duration of the walk averaged 
over the initial positions 1, . . . , N — 1 of the walker plays the role of (n) in the original 
two-walker system. 
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A standard approach for the calculation of (n) is to formulate the problem for the re- 
stricted walk between the absorbing sites in terms of a conditional first passage problem 



for an unrestricted walk on an infinite Id lattice [23]. The starting point is the Markovian 
master equation 

Pn+l(j) = \ [Pn{j - 2) + 2P n (j) + Pn(j + 2)] + [P n (j ~ 1) + P n (j + 1)] , (1) 

where P n (j) is the probability to find the walker at a given site j in the infinite lattice after 
n time steps. The first term on the right hand side of eq. (Q) is the contribution due to the 
synchronous events, by which the walker either remains at rest or it moves (here) two lattice 
sites either to the left or to the right. The second term describes jumps by one lattice site as 
a consequence of asynchronous events. The mean time to absorption (n) can now be viewed 
as a first-passage time and computed via a generating function approach |23J]. Specifically, 
(n) can be expressed in terms of the generating function for the probability of the walker 
arriving for the first time at a given r-site after a given number of steps without having been 
previously in neither of the other three r-sites. Unfortunately, the resulting expressions for 
(n) are not very transparent and their analytical dependence on N and p is not easy to 
determine. 

Another possibility is to make use of the single step probabilities appearing as coefficients 
in the right hand side of eq. (JTJ) to compute the transition probabilities between the states 
of the underlying absorbing Markov chain. In this approach, (n) can be related to the row 
sums over the elements of the fundamental matrix for the transition to the absorbing state 
or, alternatively, to its smallest eigenvalue j24[. The disadvantage of this method is that it 
requires the inversion of increasingly large matrices as N becomes large. 

Finally, the method we shall further develop here exploits the analogy of our random walk 
problem with the classical ruin problem studied by Feller [^] . Even though this approach has 
the disadvantage of being difficult to generalize to higher dimensions, it provides an elegant 
solution for the Id problem. 



3. CONNECTION WITH RUIN PROBLEM AND ANALYTICAL SOLUTION 

We first recall briefly the classical gambler's ruin problem. Consider a single walker (in 
our setting, this would correspond to the limit p=0), whose position z is viewed as the 
capital of a gambler playing against an adversary whose capital is N — z. At each time 
step, a trial is made, as a result of which the gambler wins or loses one euro. Thus, the 
gambler's winning corresponds to a nearest neighbor jump of the walker to the right, while 
losing the trial corresponds to a jump to the left. The game goes on until the gambler's 
capital is reduced to zero or increases to N (absorption of the walker at sites or N). One 
is interested in the mean duration of the game (n) z when the gambler starts with a given 
capital z. This quantity can be shown to be finite as long as < oo and obeys the following 
difference equation p: 

(n) z = ^(n) z +i + -<n)*_! + l, < z < N (2) 
with the boundary conditions 

(n)o = and (n)jy = 0. (3) 
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Eq. (J2J states that the walker has no memory of where it has been at earlier time steps; 
if it is initially at site z, it will either jump to site z + 1 or to site z — 1 with probability 
1/2. Once at any of these sites, the walker will continue his walk without remembering its 
previous position z. It is as though the walker started a new walk from z + 1 or z — 1 with 
equal probability, except that the expected value of the mean time to absorption must be 
increased by one unit j25|. The boundary conditions reflect the fact that, if the walker 
is initially placed at a r-site, it is immediately absorbed. In the original two-walker system, 
this is equivalent to placing both walkers at the same site. 

Eqs. ©-© can be solved by standard methods, e.g. by writing the general solution as 
the sum of the general solution of the corresponding homogeneous equation plus a particular 
solution [26]. One obtains in this way 



(n) z = z(N-z), (4) 

and 



z- 



The opposite limit of the above (p = 1), corresponding in our setting to a purely syn- 
chronous motion of two walkers, is somewhat less standard. In this case, only jumps by zero 
or two lattice sites may occur. Depending on its initial position, the walker may land on 
any of the four r-sites depicted in Fig. |3] Therefore, two additional boundary conditions for 
the absorbing sites —1 and N + 1 are needed. 

In the gambler's jargon, a jump by two sites means doubling the stake of each trial, 
i.e., the player wins or loses two euros with probability 1/4. Besides, a tie occurs with 
probability 1/2 (no jump). The game terminates when either the player or his adversary 
reaches or overreaches a capital of N euros. Now, the difference equation for the duration 
of the game is a fourth-order one: 

(n)z=\(n) s + 2 + ^(n) g + ^(n) g -. 2 + l, < z < N. (6) 

The solution of this equation requires four boundary conditions, namely 

(n)_i = 0, (n) = 0, (n)jv = and (n) N+1 = 0. (7) 

A particular solution of © is given by — \z 2 . The characteristic equation of the homogeneous 
equation has two double roots Ai 2 = ±1- Therefore, the full inhomogeneous solution reads 

(n) z = - X -z 2 + ci + c 2 z + (c 3 + c 4 z)(-l) z , (8) 

Substituting into eq. (jHJ) one obtains a system of equations for c\ to C4 whose solution yields 

(n) z = (z(N - 1 - z) + N) /2 for z odd, N odd, 

= z(N + 1 - z)/2 for z even, N odd, (9) 

= (z(N -z) + {N+ 1)) /2 for z odd, N even, 

= z(N — z)/2 for z even, N even. 
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The spatially averaged time to absorption (n) is again easily computed. One obtains 

(n) =JV(iV + l)(JV + 2)/12(JV-l) for iV even, 

= (N+ l)(JV + 3)/12 for N odd. (10) 

For notational convenience, let us rename the value of (n) obtained from eqs. (|l(Jj) as 
(n)^, and the corresponding result of Montroll for the one- walker problem (eq. (J5J)) as (n)^. 
A comparison between (n)( ) and (n)^ is shown in Fig. 0] As expected, the synchronous 
case becomes more efficient as soon as the lattice gets sufficiently large (N > 5). In the 
limit of a very large lattice, it is asymptotically twice as efficient as the asynchronous case 
[see Fig. E])]. 

It is also worth comparing the ^-distribution of the encounter time for the case p = 1 
(eq. 0), with the earlier known result for p = (eq. For odd N, the spatial profile 

in z displays some qualitative similarities in both cases (cf . Figs. and Eb) • The highest 
value of (n) z is attained at those sites z that maximize the distance d z = min(z, N — z) from 
the closest absorbing site. For p — 0, the encounter time always increases with increasing 
d z , and it becomes maximum when z = (N ± l)/2. In contrast, the behavior is no longer 
strictly monotonic for p = 1, since (n) z either increases or remains constant as d z becomes 
larger, giving rise to a "staircase" profile for high enough values of A" (Fig. [BJ}). 

For even lattices, there are more marked differences between the cases p = and p = 1. 
For p = 0, the discrete spatial profile resembles an inverted parabola like in the odd lattice 
case, the maximum now being located at z = N/2. However, if p = 1, one observes a series of 
alternating valleys and peaks in the distribution (cf Figs. [7^ and[7b). There are two subcases 
here: if A" is divisible by 4, the highest values of (n) z are attained for z = (N ± 2)/2 (Fig. 
Ek). Otherwise, the maximum value corresponds again to z = N/2 (Fig. Ujp). 

Note that for p = and even A^ the time to absorption when the walker is started at sites 
1 or A^ — 1 is smaller than in the purely synchronous case. An intuitive argument points to 
the fact that, in the former case, the walker is absorbed one out of two times after the first 
time step, while absorption takes place only one out of four times if p = 1. However, this 
argument should be taken with care, since it fails for odd N. Besides, as we shall see later 
on, the minimum of (n)i = (n)jv-i corresponds to a process with p > 0. 

We finally turn to the general case of a mixed walk. Assume that a given event is 
synchronous with probability < p < 1 and asynchronous with probability 1 — p. The 
difference equation for (n) z now reads: 



(n) z = |(n) 2+2 + ^l(n) z+1 + |(n) a + —^■(n) g - 1 + ^(n},_ 2 + 1, < z < N. (11) 

The boundary conditions are again given by (JJJ). A particular solution of eq. (fTT|) is given 
by — z 2 /(l +p). The roots of the underlying characteristic equation are 

Al , 2 = l, A» = ± ^HZ . (12) 

P 

The full solution reads 

(n) z = -z 2 /(l +p) + c 1 + c 2 z + c 3 A3 + c 4 X z 3 , (13) 
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where the constants c\ — c 4 are now given by the linear system 



h Ci - c 2 + A 3 1 c A + A 4 x c 4 = 0, (14a) 



1 +p 

ci + c 3 + c 4 = 0, (14b) 
iV 2 

+ ci + Nc 2 + \%c 3 + Af c 4 = 0, (14c) 



1 +p 
(N+l) 2 

l+P 



+ a + (N + 1) c 2 + A^ 1 c 3 + Af +1 c 4 = 0. (14d) 



We prefer to omit the rather lengthy analytic form of the coefficients C1-C4. From the above 
equations, (n) z and (n) can be explicitly computed. For (n) one has an expression of the 
form 

{U} ~ ~ 6(1 +p) +C ^ N >P> + 2 

M^p) [Asfr)]"" 1 + c 4 (iV,p) [A 4 (p)f - 1 ) . (15) 

Tabled displays the analytical expressions of (n) as a function of p for increasing values of iV. 
These are rational functions of p whose complexity increases with N. For N > 4 and small 
p, one can obtain a first-order approximation to the exact solution by Taylor-expanding eq. 

^ + ^p + o(p 2 ) 



(n) = (n) 



(o) 



(16) 



Figs. |St and |Hb contain several numerical plots of (n), as a function of p for different 
values of N. For N = 2, the purely asynchronous case is more effective than the purely 
synchronous case 1 . As synchronicity is turned on, a monotonic increase of (n) as a function 
of p is observed. For N = 3, (n) does not depend on p. At each individual time step, the 
probability of reaction is always 1/2, regardless of whether both walkers hop or only one of 
them. In a sense, the cases N = 2 and N = 3 are non-generic, since they only involve a 
single symmetry- distinct initial condition. For iV > 4, a new parity effect appears: if iV is 
even, the most effective process is observed for an intermediate value of p associated to a 
minimum of the function (n) in the physically acceptable p-interval [0,1]; in contrast, for 
odd lattices, (n) is a monotonically decreasing function of p, i.e., the most effective process is 
always the purely synchronous one. For N = A, any intermediate value of p makes hopping 
more effective than in both limiting cases, and a minimum of (n) is obtained for p m i n = 2/3. 
For higher, even values of N, the minimum is rapidly shifted to the right (p m in ~ 0.86 for 
iV = 6) and the p-interval for which processes are more efficient than the purely synchronous 
case shrinks dramatically. For large N, p m i n gets arbitrarily close to 1 (see Table E]). 

A series of Monte Carlo simulations for the periodic two-walker system has been carried 
out to confirm our analytic results based on the one-walker description. For two different 
lattice sizes, namely N = 7 and iV = 8, we have performed a series of statistical runs to 
compute (n), each run thereby comprising a whole set of symmetric-distinct nonreactive 



Obviously, the purely asynchronous case has a maximum efficiency in this case, since the walker is always 
trapped after the first step. 
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configurations. Due to the large variability (n 2 ) — (n) 2 characteristic of first-passage prob- 
lems, a relatively high number of runs was needed (10 6 ) to obtain an accurate value for (n). 
The (n)-values obtained from simulations (with an accuracy of three significant digits) are 
listed in Tables El and |U for the cases N = 7 and N = 8, respectively The agreement with 
the analytical predictions is good, the maximum observed deviation is off the theoretical 
value by about 1% only. 

In order to obtain additional insight in the even-odd transition mechanism, we have 
studied qualitatively (n) z as a function of p and iV for each single initial position z of the 
walker. Let us characterize each site z by its distance d z to the closest r-site. For not too 
large values of p, the behavior of the encounter time is roughly the same for all z values 
regardless of the parity of N, i.e., a decrease of the encounter time is observed (Figs. and 
Et>). However, the qualitative p-dependence of (n) z in the large p limit becomes different 
for initial positions with even or odd values of d z : for even values of N and initial positions 
with odd d z , (n) z begins to increase sharply, as a result of which a minimum of the curve is 
observed (cf Fig. EK)- Even though the contribution to the global efficiency (n) arising from 
the {N/2) — 1 sites with even d z decreases strongly in this regime, this effect is overcome by 
the increase of the contribution yielded by the N/2 sites with odd d z , thus giving rise to a 
net increase in (n). In contrast, the sensitivity to the initial condition is less systematic and 
less important for odd values of N (Fig. again, a monotonic decrease is observed for 
even d z . Even though minima are still observed for d z = 1, for all other odd values of d z , 
they are either absent 2 or very flat (see curve for d z = 3 in Fig. Eb)- 



4. MONTE CARLO RESULTS IN TWO AND THREE DIMENSIONS 

To complement the analytic Id results for the dependence of (n) on the value of the 
parameter p, we have also investigated how (n) depends on p for dimensions 2d and 3d using 
Monte Carlo simulations. Figures and show the dependence of (n) with respect to p. 
As shown above for Id, there is distinctively different behavior based on whether N is even or 
odd. In higher dimensions this even-odd effect becomes much more pronounced. It is clear 
that when N is an even number, the maximal efficiency is attained at some intermediate 
value between and 1. Also, one notices that for the 6x6 square lattice and for all 3d even 
lattices which we studied (n)^ < (n)^. This is a surprising result which we do not see in 
the Id case for lattices with N > 4. This suggests the existence of a crossover effect in the 
large size limit when switching from 2d to 3d lattices, implying that in the former case the 
two simultaneously moving walkers are more effective than two asynchronous walkers, while 
in 3d the opposite holds. As in the Id case, the value of p min in 2d and 3d tends toward 1 
in the limit of large lattice size, but it is interesting to note that in the largest 3d lattice 
which we studied (N = 1000), the difference between p = 0.999 and p = 1 is ~ 600, which 
is about 30% of (n)^ in that case. This shows that for even lattices, a minute amount of 
asynchronicity allows for a much greater efficiency Similar arguments to those given above 
for Id exist for in 2d and 3d when attempting to determine why such even-odd behavior 
arises, and further analysis of this striking behavior is given in Ref. |27| . 



2 This is e.g. the case for N = 7 and d z = 3 (not shown here). 
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5. COMPARISON WITH CONTINUUM APPROXIMATION 



It is instructive to compare the above results in Id with the continuum approximation 
valid for large N. To do so, consider the one-walker system with absorbing sites and a fixed 
lattice length L = NAx, where Ax is the intersite distance (lattice constant). According 
to this definition, L is the distance between the inmost r-sites. The walker's distance to 
site is x = z Ax. We perform the continuum limit by letting Ax and At simultaneously 
go to zero under the additional requirement that the diffusive combination (Ax) 2 / At tend 
to a finite constant. Since L is fixed, this implies that one lets the number of sites N go 
to infinity while the spatial and temporal resolutions Ax and At are scaled as 1/N and 
1/iV 2 , respectively. Let us next replace (n) z by a function (n) x varying smoothly in the 
space interval [0, L\. The mean elapsed time (t) x to absorption will then simply be the mean 
number of steps (n) x times the time unit At. From eq. (fTTj) . we have: 

- [{t)x+2Ax ~ 2(t) x + (t) x -2Ax} + [(t) x +Ax - 2(t) x + (t) x ^ Ax \ +At = (17) 

with the boundary conditions 3 

(*)-A* = (t)x = (t) L = (t) L+Ax = 0. (18) 

We now divide eq. (|17|) by At and expand the expressions in the brackets in Ax. Taking 
the diffusive limit in the resulting equation yields 



0^ = -i, (19) 



where the relative diffusion coefficient D is given by 

D = (l+p) lim ^L = (i+p)£)(°). (20) 
v ry As,At— o At 

In the rightmost equation, is the value of the diffusion coefficient for p = 0. For p > 0, 
is increased by the prefactor 1 + p, i.e., the variance of the single-step probabilities of 
the random walk. In this limit, the four boundary conditions (fTHj) coalesce into two distinct 
ones, namely (t) = and (t)i = 0. The solution of (fT^j) which fulfils these boundary 
conditions is 

W. = X ~^- (21) 
The spatially averaged reaction time (t) is obtained by integrating over x: 

<«> = ir< t >" fa =i^- (22) 

As expected, (t) is proportional to the squared lattice length and inversely proportional to 
the relative diffusion coefficient D. For the special cases p = and p = 1, this result is 



In the case p = 0, one only has two boundary conditions, namely (t)o = (t) l = 0. 
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recovered by directly taking the diffusive limit in the discrete solutions (J3J) and (|1U|). The 
relation (}2"2*j) leads to the asymptotic law 

/+\(0) / n \(o) 

%- = lim ^V = l+p, (23) 

(i) v-^oo (n) 

where (t) (0) = L 2 /(12D^) and (n) is given by eq. ([T5|). Eq. (H3J) shows that in the 
continuum limit, the efficiency of the reaction increases linearly with p. The continuum 
approximation applies for sufficiently large N, namely when the typical displacement of the 
walker Al = ^1 + p Ax at each time step is small compared to the lattice length NAx. If 
this condition is not fulfilled, the approximation gets significantly worse at large values of p 
(cf Fig. H2J). 

A generalization of (|19|) for higher order moments {P) x can be obtained by writing down 
the difference equations for the discrete quantities {n?) z and taking the diffusive limit thereof. 
One then gets the coupled set of equations 

d 2 (t j+1 ) 

D ~^ h = -^ + l )^ ( 24 ) 

Eqs. (124}) . not to be further dealt with here, are well known from the theory of first-passage 
problems, where they are usually obtained from the adjoint Fokker-Planck equation for the 
underlying diffusion process j2^,|2^. Again, deviations from the dynamics dictated by 
are expected for small lattices. 



6. CONCLUSIONS AND OUTLOOK 

We have seen that the Id problem of computing the mean reaction time between two 
diffusing co-reactants can be reduced to a trapping problem for a single walker. The latter 
can be viewed as a generalized ruin problem, the duration of the game plays thereby the 
role of the mean time to absorption. 

In the diffusive limit, equivalent to the limit N —>■ oo if the lattice length L is held fixed, 
the reaction efficiency increases linearly with p, but important deviations are observed for not 
too large values of N. Beyond the crossover size N = 4, a new parity effect is observed. For 
odd values of N, the reaction time still increases monotonically (but no longer linearly) with 
p, while for even values of N, the efficiency is optimized for an intermediate value p m % n < !■ 
In higher dimensions, this parity effect is even more pronounced, i.e. for even lattices there 
is a drastic increase in efficiency when a tiny amount of synchronicity is introduced. In 
contrast to the Id case, the effect is enhanced with increasing N. 

Let us briefly comment on the Id results from the perspective of the ruin problem. Assume 
that one of the gamblers is successively given an starting capital of 1, 2, . . . , N — 1 euros at 
each round, while his adversary gets N — 1,N — 2,...,1 euros. Let us further suppose that 
the gamblers can choose between two kinds of trials: the stake for the first one is one euro 
and there is no tie. In the second trial, either nobody wins or one of the gamblers wins two 
euros. According to our results, for an odd capital N, the gamblers minimize the average 
playing time if they always make a two-euro bet. However, if N is even, they should make 
a small amount of one-euro bets in order to finish the round as soon as possible. 

Our work can be extended in many different ways. Perhaps the most straightforward 
one is the detailed characterization of the whole distribution P(n) in terms of p and N. 
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A generalization of results such as equation (|23|) to higher dimensional integral and fractal 
lattices is also of interest, since it may further clarify the role of dimensionality and the 
lattice coordination number. According to our results, in one-dimension two synchronously 
moving walkers are asymptotically twice as efficient as when they hop one after the other. In 
a 2d square planar lattice, Kozak et al. have shown that the purely synchronous case is v2 
times more efficient than the purely asynchronous one in the large N limit for odd lattices 
[lfl| . The question is whether or not the relative efficiency of both processes in lattices with 
fractal dimension 1 < df < 2 lies between 2 and a/2. Preliminary calculations on a Sierpinski 
gasket (with fractal dimension df — 1.585) seem to indicate an asymptotic relative efficiency 
higher than 2 in this case, despite the fact that the lattice has (up to the three vertex sites) 
the same coordination number as a 2d square planar lattice [l?}. The reason for this may be 
the important role played by the specific form of the lattice boundaries, even in the limit of a 
large lattice. This may motivate the study of boundary conditions other than periodic ones 
for the two-walker system. However, the analytical treatment of this case is considerably 
harder, at least in the framework of the method of difference equations, since the lattice is 
no longer translationally invariant. 

As a further extension of our work, one can also consider more complex reactive schemes 
[29! ] involving more than two walkers to study the combined effect of synchronicity and 
many-body effects. 
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Figure Captions 

Figure 1: a) two- walker system on a seven-site periodic lattice (walkers represented by black 
circles), b) Equivalent one- walker system with trap. For convenience, the walker labels A 
and B in Fig. ^ have been left out. The arrows in Fig. ^ indicate that both walkers 
perform a synchronous step. In Fig. lb this corresponds to a two-site jump of the walker. 

Figure 2: Lattice transformation for the one-walker system displayed in Fig. ^p. 

Figure 3: Replacement of the traps T at each end site of the transformed lattice by two 
absorbing sites r. 

Figure 4: Mean encounter time as a function of the lattice size for the purely synchronous 
case (n) = (n)^ 1 ' (circles) and the purely asynchronous case (n) = (n)^ (crosses). 

Figure 5: Ratio (n)^ / (up 1 * as a function of N. Note that the value of the ratio is the same 
for two consecutive odd and even values of N. The inset displays the behaviour for small 
N. 

Figure 6: ^-distribution of the encounter time in the cases p = (dashed lines) and p = 1 
(continuous lines) for a) iV = 3, 5, 7 and b) N = 9, 11. For N = 3 both cases display the 
same flat distribution. 

Figure 7: z-distribution of the encounter time for p = (dashed lines) and p = 1 (continuous 
lines) for a) N = 6, 10 and b) iV = 4, 8. 

Figure 8: Mean encounter time as a function of p for a) N = 2, ... ,5 and b) N = 6, . . . ,9. 

Figure 9: a) p-dependence of (n) z for all possible values of d z on a lattice with a) 10 sites 
b) 9 sites. 

Figure 10: Mean encounter time as a function of p for two walkers on a 2d square planar 
lattice with periodic boundary conditions. 

Figure 11: Mean encounter time as a function of p for two walkers on a 3d cubic lattice with 
periodic boundary conditions. 

Figure 12: Ratio (n)^/(n) as a function of p. 
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5 4(2p- 5)/(p 2 -4) 

6 (28/5) (p 2 - 10p + 10) /V - 4p 2 - 4p + 8) 

7 (4/3) (p 2 + 8p- 14) /(p 2 - 2) 

8 (12/7)(13p 3 + 6p 2 - 126p + 112)/((p - 2)(p 3 + 6p 2 - 8)) 

9 10(2p 3 - 5p 2 - 16p + 24)/((p 2 + 2p - 4)(p 2 - 2p - 4)) 

10 (22/9) (7p 4 - 76p 3 + 16p 2 + 288p - 240)/(p 5 - 6p 4 - Up 3 + 32p 2 + 16p - 32) 



TABLE 1: 


Analytic expressions for (n). 


iV 




2 


0.0000 


4 


0.6667 


6 


0.8596 


8 


0.9204 


10 


0.9483 


12 


0.9636 


14 


0.9729 


16 


0.9791 



TABLE 2: Values of Pmin with 4-digit accuracy. 
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0.2 ~ 8.410 8.408163 

0.5 7.430 7.428571 

0.8 6.826 6.823529 

1 6.667 6.666667 



TABLE 3: Theoretical vs. simulation value of (n) for N = 7. 



p ^ n ^sim / n \theo 

0.2 10.709 10.706177 

0.5 9.345 9.344538 

0.8 8.495 8.496241 

1 8.572 8.571429 



TABLE 4: Theoretical vs. simulation value of (n) for iV = 8. 
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